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Abstract 

In this paper we discuss the production of a heavy scalar MSSM Higgs boson H 
and its subsequent decays into pairs of electroweak gauge bosons WW and ZZ. 
We perform a scan over the relevant MSSM parameters, using constraints from 
direct Higgs searches and several low-energy observables. We then compare the 
possible size of the pp — > H — > WW, ZZ cross sections with corresponding Stan- 
dard Model cross sections. We also include the full MSSM vertex corrections to 
the H — > WW, ZZ decay and combine them with the Higgs propagator correc- 
tions, paying special attention to the IR-divergent contributions. We find that 
the vertex corrections can be as large as —30% in MSSM parameter space regions 
which are currently probed by Higgs searches at the LHC. Once the sensitivity of 
these searches reaches two percent of the SM signal strength the vertex correc- 
tions can be numerically as important as the leading order and Higgs self-energy 
corrections and have to be considered when setting limits on MSSM parameters. 
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1 Introduction 



The recent discovery of a 126 GeV resonance decaying into photons and (off-shell) Z 
bosons at the LHC [1, 2] opens a new era in particle physics. The next important task for 
both theorists and experimentalists is to determine the exact nature of that resonance. 
Currently the measured signals are in statistical agreement with the expectations from 
a Standard Model (SM) Higgs boson. However, the experimental sensitivity is not yet 
sufficient to rule out an extended Higgs sector, especially if the (tree-level) couplings of 
the additional Higgs bosons to electroweak gauge bosons are suppressed. The discovery 
of additional scalar resonances would give us important clues about the exact mechanism 
of electroweak symmetry breaking. 

In the minimal supersymmetric SM (MSSM) the tree-level couplings of the neutral Higgs 
bosons to weak gauge bosons are determined by the Higgs mass scale (either Ma, the 
mass of the pseudoscalar Higgs boson, or M H ±, the mass of the charged Higgs boson) 
and tan/3, the ratio of the vacuum expectation values of the Higgs doublets. At leading 
order the masses of the Higgs bosons are also determined by these two parameters. For 
Ma ^> Mz, the so-called decoupling limit [3, 4], the heavy scalar Higgs boson H and 
the pseudoscalar A are almost degenerate and their (effective) couplings to W and Z 
bosons are strongly suppressed. This makes the search for heavy MSSM Higgs bosons 
more difficult than the search for a Standard Model Higgs boson with similar mass. 
However, it is well-known that the masses and couplings of MSSM Higgs bosons receive 
large corrections at higher orders in perturbation theory [5-8]. Also, the production 
rates for Higgs bosons are modified in the MSSM, especially in the gg — >■ H, A and 
bb — )> H, A production modes. 

In [9] the production and decays of a pseudoscalar Higgs into electroweak gauge bosons 
were discussed in a number of different models, including the MSSM. In this paper we 
answer the question of how large the LHC signal cross sections for pp — >■ H — >■ WW, ZZ 
can become in the MSSM when higher order corrections to both the production and 
decay processes are taken into account. For this purpose we perform a scan over the 
relevant MSSM parameters, using experimental constraints from several low-energy ob- 
servables and direct Higgs searches at LEP, Tevatron and LHC. We do not assume a 
specific SUSY breaking scenario, but scan directly over the soft SUSY breaking param- 
eters at the electroweak scale. For this scan we make extensive use of the public codes 
HiggsBounds 3.8.0 [10, 11] and FeynHiggs 2.7.4 [12-18]. 

The Higgs-gauge-boson couplings are implemented in FeynHiggs in the improved Born- 
approximation, i.e. taking into account higher order corrections from Higgs self-energies 
but no genuine vertex corrections. The MSSM vertex corrections for both the WW and 
ZZ final state were calculated in [19], although for the WW final state only fermion and 
sfermion contributions were considered. For our analysis we performed an independent 
calculation of all one-loop vertex corrections and found agreement with [19]. We then 
extended the analysis of the H — > WW case to the complete MSSM corrections, in- 
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eluding the IR divergent contributions and the corresponding real emission graphs. Our 
scan shows that the vertex corrections typically lie between —10% and —30% in MSSM 
parameter space regions where the H — > WW, ZZ channels should still be observable 
at the LHC. 

The case of off-shell decays of the light MSSM Higgs-boson h was discussed in [20], where 
a calculation of the process h — > W*W*, Z*Z* — > 4f (four fermions) was presented. In 
this paper we examine the possibility of calculating the single off-shell processes H — > 
WW* — > Wff and H — > ZZ* — > Wff process in an effective coupling approximation, 
i.e. by re-scaling the corresponding SM decay rates. Such an approximation can be 
useful in parameter scans or fits, where off-shell decays of the heavy MSSM Higgs boson 
may be of interest, but a numerical integration of the full four-particle phase space is not 
feasible. We discuss the quality of the approximation and address the issue of infrared 
divergences in this approach. 

In Section 2 we introduce our notation and explain the combination of the vertex cor- 
rections with the self-energy corrections calculated by FeynHiggs. In Section 3 we give 
the details of the parameter scan and discuss the experimental constraints that were 
used in it. The numerical results of the scan and the quality of the effective coupling 
approximation are discussed in Section 4. Our conclusions are given in Section 5. 

2 Details of the calculation 

2.1 Notation and conventions 

In the MSSM, the Higgs sector contains two scalar doublets, which give five physi- 
cal Higgs bosons. At lowest order, the Higgs sector is CP-conserving, containing two 
charged Higgs bosons, H ± , two neutral CP-even Higgs bosons, h and H, and the CP-odd 
Higgs A. Two independent parameters characterise the Higgs sector, normally taken as 
Ma and tan/3, where tan/3 is the ratio of the vacuum expectation values of the Higgs 
doublets. Higher order corrections lead to large corrections to the Higgs masses and 
mixing angle a, and can induce CP- violation and mixing between the three neutral 
Higgs bosons h, H and A [21-23] if complex SUSY-breaking parameters are allowed. 1 

2.2 Higgs propagator corrections 

Higgs propagator corrections can be extremely important numerically, especially in the 
non-decoupling regions of the SUSY parameter space, and are in addition needed in 
order to ensure correct on-shell properties of S'-matrix elements involving external Higgs 

1 In the case of CP-violation, it is usual to take M H ± as input parameter instead of Ma because, in 
the CP-violating case the pseudoscalar Higgs boson A mixes with the CP-even neutral Higgs bosons. 
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bosons - i.e. unit residue and vanishing mixing between different Higgs bosons on mass 
shell. These corrections can be included by using finite wave function normalisation 
factors. In the case where these factors are applied to a tree level decay amplitude we 
speak of an improved Born approximation. In the following, quantities computed in 
this approximation are denoted with a subscript 'imp.B'. An amplitude Af/.imp.B in the 
improved Born approximation with an external Higgs boson H can receive corrections 
from three tree- level amplitudes Ah,tree, Aff,tree and v4.A,tree involving the three neutral 
Higgs states: 

Af^imp.B = Zjjfc.A^tree + %HH Ah, tree + ^HaAa^tcc (1) 

The matrix Z has been defined in Ref. [15, 24] and is non-unitary. When no CP- 
violation is present mixing occurs only between the CP-even states, but when complex 
parameters are allowed mixing between all three neutral states needs to be considered. 2 
The program FeynHiggs 2.7.4 [12-18] has been used to calculate both the corrected 
Higgs boson masses and the wave function normalisation Z-factors. FeynHiggs includes 
the complete one-loop corrections as well as the dominant two-loop contributions in the 
MSSM with real and complex parameters. 

Since the Higgs propagator corrections are universal, they can in principle be applied 
to the loop diagrams as well as the tree-level diagrams. Denoting the one-loop vertex 
corrections to the decay amplitudes of the tree-level mass eigenstates h, H and A as 
A Ah, A Ah and AAa, respectively, we define the improved vertex corrections for the 
physical mass eigenstate as 

A Amp = Z Hh AA h + ZhhAAh + ZhaAAa ■ (2) 

When computing interferences between the tree-level and one-loop vertex diagrams, 
improved versions can be used for neither, the tree-level or both of the factors. This 
provides an easy method of including (potentially large) higher-order corrections in our 
calculations. In the CP-violating case, applying the propagator corrections at loop level 
could give rise to interesting effects as it allows the CP-odd Higgs boson, A (which of 
course does not couple to the gauge bosons at tree level), to be taken into account. 

In any case, applying the Higgs propagator corrections means that we are mixing pertur- 
bative orders and could potentially miss cancellations found at higher orders. However, 
estimations of the uncertainties from unknown higher order corrections (see [7, 8, 25]) 
indicate that the Z-factors do indeed give rise to a leading contribution which is not 
expected to be numerically compensated by the remaining 2-loop pieces. Since the ef- 
fect of applying the Higgs propagator corrections at loop level is significant (as we shall 
show), we choose to follow this method. 

2 For the H — > VV decays considered in this paper, _4^ co is of course zero. 
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2.3 Comparing the SM and MSSM 



In this paper we are interested in MSSM scenarios that lead to relatively large pp — > 
H — > VV signals (VV = WW, ZZ). Hence, we define the ratios 

_ a™ SSM • BR(H -> VV) MSSM _ a™ SSM rf 4 



]•(// > yy) mssm 

t(h ->• vvy 



\ — ' y v ) /o\ 



where a// denotes the LHC production cross section for the Higgs H, BR(H — > VV) 
and T(H — > VV) are the branching ratio and partial decay width into vector bosons V 
(V = W, Z) and Y H is the total decay width of the Higgs boson H. The superscript 
'MSSM' indicates that the corresponding quantity is evaluated in the MSSM, with 
H being the heavy scalar MSSM Higgs boson. The superscript 'SM' means that the 
quantity is evaluated in the SM, with H being a SM Higgs boson with the same mass 
as the heavy MSSM Higgs boson. At leading order the ratios Ryv and pvv are the 
same for V = W and V — Z, since the ratio between the Standard Model and MSSM 
couplings is the same for both HWW and HZZ. From the definitions of Eq. (3) it is 
obvious that the pp — > H — >■ VV cross sections and H — > VV partial widths within the 
MSSM can be obtained by scaling the corresponding SM quantities with R vv or p V v- 

If the Higgs mass Mh is below the VV threshold the Higgs boson H may still decay 
into Vff (with / and /' being light fermions) via an off-shell vector boson V*. If 
new-physics contributions to the Vff vertex and non-factorisable contributions are 
neglected, the corresponding ratios of partial widths or cross sections times branching 
ratios do not depend on the fermions / and /'. For off-shell decays we therefore define 
Rvff and pvff as ratios of differential cross sections and partial widths: 



PVff'i M ff') 



dT(H -)■ VV* -)■ Vff) MbbM /dM ff 
dT(H VV* Vff') SM /dM ff , 



0.MSSM pSM 

R VJf ,{M fr ) = -*.^^.p vjr {M Jf .) , (4) 
h a 

where Mff denotes the invariant mass of the //' pair. Differential pp — >■ H —¥ Vff 
cross sections time branching ratios and differential H — > Vff partial widths within 
the MSSM may thus be obtained by scaling the corresponding SM quantities with Rvff 
and pvff - Usually, pvff is only weakly dependent on Mff. We may then approximate 
pvff(Mfft) as a constant, 

Pvff'(M f y) w Pvff(M H - M v ) = p vf f , (5) 

and calculate integrated MSSM cross sections and partial widths by scaling correspond- 
ing SM quantities with the appropriate factors. This approximation is what we call the 
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effective coupling approximation, since higher order corrections to the HVV vertex have 
been absorbed into an effective coupling constant. 

The principle behind this effective coupling approximation is the same as that used 
by the Higgs Cross Section Working Group when working in the MSSM. In order to 
include all known higher order corrections (some of which are known only in the SM, 
not the MSSM), the Working Group takes SM 'building blocks' and dresses them with 
the appropriate MSSM coupling factors, as described in [26]. 



2.4 Higher-order corrections and form factors 

We can incorporate the corrections to the HVV vertex by calculating an effective HVV 
coupling resulting from the loop and counterterm diagrams. The structure of this cou- 
pling for on-shell particles is [27-29] 

T^( qi , q 2 ) = A( Ql , qfigT + B(q u + C(q u q 2 )e^q lp q 2r7 . (6) 

Here, qi and q 2 are the momenta of the electroweak gauge bosons, and A, B and C 
are Lorentz invariant form factors. For off-shell particles, the coupling can have a more 
complicated structure, but if the gauge bosons decay into massless fermions the only 
relevant form factors are A, B and C. At tree level, only the formfactor A has a non-zero 
value in both the SM and the MSSM: 

SM _ ieM w SM _ ieM w 

A »sMww- sin6w , a HsmZZ - s]ndwCOs2dw in 

A^ M = A™ MVV sm(P-*) , A™™ = A™ mVV cos((3 - a) ,VV = WW,ZZ 

(8) 



At lowest order the MSSM formfactor A representing the coupling of the light CP-even 
Higgs boson differs from the SM value of A by a factor of sin (/? — «), which tends to 1 
in the decoupling regime, i.e. for Ma ^> Mz. Higher order diagrams, however, lead 
to different contributions to A in the Standard Model and MSSM, and can result in 
non-zero values for B and C. 

For the calculation of the form factors we employ a mixed renormalisation scheme where 
the electroweak sector is renormalised on-shell [30] , while the Higgs sector is renormalised 
using a hybrid scheme where the Higgs fields are renormalised in the DR scheme and 
Ma is renormalised on-shell, as described in [15]. We parameterise our results in terms 
of a{Mz) and calculate the charge renormalisation constant accordingly - i.e. 

5Z e -> 5Z e - ^Aa (9) 

For Higgs bosons inside loops we use the physical masses and the unitary Higgs mixing 
matrix calculated by FeynHiggs, as described in [15]. 
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2.4.1 Infrared divergences in H — > WW 



In the tensorial structure given in Eq. (6), it is only the form factor A that is IR- 
divergent due to photon exchange. To render transition probabilities finite one must 
sum over all energy-degenerate final and initial states [31, 32]. In practice this means 
that one has to include contributions involving real radiation of a photon in order to 
obtain infrared-finite observables. In the SM the analytic expression for the on-shell 
real correction to the partial width H S m — > W + W~ reads 



T 3 (H SM ^WW)=T - 

7T 



-F 1 (/3 )log^ + 2(2-log[2]) + i±^ ^ 2 '" 



Li 2 —^- (10) 



o 

] j 2 >'<> T: I ""A) I; 1 - A. , 1 , ri , 2ll I - >0 



27^- + Li 2 ^^ - Li 2 - — '— + - log[l - B'} log 
B Q + 1 2 2 2 B[ Hoi l-f3 

M% a a ( / a n a x 1 + Bn K 4, 



with 



and Mh sm the mass of the SM Higgs boson. The function Fi(Bq) is given by 

W) = i±^logi±A-2 , (12) 

Po J- — Po 

and r denotes the lowest order SM decay width, given by 

When T 3 {H — > WW) is added to the correction to the partial width from loop diagrams, 
r virt (ff — > WW), the virtual corrections cancel the IR divergences in the real corrections 
order-by-order, as required by the Bloch-Nordsieck theorem [33]. 

To improve the accuracy of our results in the MSSM we would like to apply the Higgs 
propagator Z-factors to the vertex diagrams as well as the tree diagrams, as described 
in Section 2.2. However, as discussed previously, we are then mixing different orders 
of perturbation theory and the cancellation of IR divergences is no longer guaranteed. 
Once the Higgs propagator corrections have been applied according to Eq. (1), the 
leading order SM and MSSM couplings are related through the following equation: 



A imp . B = (sin(/3 - a)Z Hh + cos(/3 - a)Z HH ) A 



SM 
tree 



"MSSM 



AH ■ (14) 



Here and for the rest of this section, the superscript 'SM' indicates that the correspond- 
ing quantity is evaluated in the SM with the SM Higgs boson mass set to the mass 
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(a) (b) 



(b) IR-divergent contributions to the renormal- (c) Real corrections to the HWW vertex. 

isation constants contained within the HWW 

counterterm. 

Figure 1: IR divergent diagrams contributing to the corrections to the H — > WW partial 
width. 

of the heavy MSSM Higgs boson H . Symbols without superscripts refer to the MSSM 
unless stated otherwise. When propagator-type corrections are applied to the loop and 
real radiation diagrams, the IR divergent NLO diagrams (shown in Figure 1) are also 
modified. Counterterm diagrams (involving renormalisation constants with IR diver- 
gent contributions, as shown in Fig. 1(b)) and the real radiation diagrams (Fig. 1(c)) 
are modified by the factor -F MSSM , as are most of the loop diagrams of Fig. 1(a). The 
only exception is diagram (b) of Fig. 1(a), which involves the coupling between a neutral 
Higgs and a pair of charged Goldstone bosons: the SM and MSSM couplings are related 
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by 3 : 

Ahgg = (sin(/3 - a)m 2 htrcc Z Hh + cos(/3 - a)m 2 Htrcc Z HH ) A S ^ G , (15) 

Hsm 

- i-e. Ahqq ^ F MSSM A^ G . Since m^ tree 7^ m# 

,tree 7^ "^i?sM ' ^ e diagrams involving 
the coupling between a neutral Higgs and a pair of charged Goldstone bosons in the SM 
and MSSM are not related by the same factor as the other IR divergent diagrams (or 
the real correction diagrams) and the IR divergences therefore do not cancel between 
the real and virtual contributions when Higgs propagator type corrections are applied 
at loop level as well as at tree level. 

By keeping the corrections strictly at the one-loop level this problem can of course 
be avoided, resulting in an IR finite result, just as in the SM. In this approach, an 
improved Born approximation is used for the "leading order" form factor Amp.B ~~ i- e - 
the propagator-type corrections are applied to the tree level form factors in the following 
manner: 

I^NLOI 2 = |Amp. B | 2 + 2 Re K*rcc A ^] + ^real | ^-tree I 2 , (16) 

where A tree is the tree level MSSM formfactor without propagator factors, AA is the 
correction to the form factor A arising from the virtual MSSM corrections (Z-factors are 
not applied to the loop diagrams), and <5 rca i | A tree \ 2 is the correction to the form factor 
resulting from the real radiation 4 . While this approach does avoid the problem with 
IR divergences, it has a drawback because it misses the potentially large corrections 
arising from Higgs mixing at the loop level (and only gives a correction for the CP-even 
Higgs boson decays). Several alternative approaches have been investigated, to allow 
the Higgs propagator type corrections to be included at loop level as well as at leading 
order while preserving an IR-finite result. 

• Option 1: Strictly speaking, the IR divergences are a higher-order effect - they 
occur only because we are mixing orders by applying the Z-factors at the one- 
loop level. The IR divergent terms can therefore be calculated analytically and 
subtracted "by hand" : 

|^4imp.NLo| 2 = |^imp.B| 2 + 2 Re [A* mfxB AAj mp ] + 5 r cal | ^imp.B | 2 — ^sub l^imp.B 1 2 j (17) 

where AA; mp is the contribution to the form factor from vertex corrections to 
the decay of the Higgs boson (with Z-factors applied to the loop diagrams), and 
4ub|Amp.B| 2 is the analytically-calculated subtraction term used to ensure that 
the squared form factor is IR-finite. 

• Option 2: A second approach is to treat the 'problematic' loop diagram (and 
the corresponding part of the real radiation) strictly at one-loop level - no Z- 
factors are applied to this part of the correction - whilst applying Higgs propagator 

3 Note that the a used here is the tree-level value, rather than the so-called 'effective' a often used 
to account for Higgs mixing. 

4 Throughout, a capital A implies an absolute correction and small S a relative correction. 
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corrections to all other higher order diagrams. In this case, 



|^4imp.NLo| 2 — |^imp.B| 2 + 2 Re [A* mp B AA; mp ] + <5 rca i | Ai mp 3 1 2 

+ 2 Re [A t * rcc A^ oldstonc ] + 5f^ stonc \A tree \ 2 , (18) 

where A74 goldstone is the contribution to the formfactor from the virtual correction 
shown in Fig. l(a)(b) 5 , 5f c °af stonc is the corresponding part of the real radiation 
and AA' imp is the contribution from all virtual corrections other than the loop 
diagram containing the HG + G~ coupling, with Z-factors applied. 

• Option 3: The origin of the remaining IR divergences is a mismatch between 
the HGG coupling and the Higgs mass eigenstates. By applying Z-factors, the 
latter are determined at the one-loop level while the HGG coupling is taken at 
tree-level. In a complete two-loop calculation this mismatch would be cured by 
one-loop corrections to the HGG vertex. We therefore work with an effective HGG 
coupling, whose value is completely fixed by the requirement that IR divergences 
cancel at the one-loop level: 

A »gg = ~ 9M J* a C^Hh sin(/3 - a) + Z HH cos(/3 - a)) 

ZM\y Sill u\y 

= F M ssmAh sm gg ■ (19) 
With this effective coupling, the NLO HVV form factor is 

l^imp. NLO | — |Amp.B | 2 + 2 Re [A* mp . B AAmp] + 5rcal| Amp.Bl 2 , (20) 

where AA imp denotes the one-loop vertex corrections with Z-factors applied and 
using the effective HGG coupling from (19). 

The above options are compared in Figure 2, which shows the relative correction to 
the formfactor \A\ 2 (as given by Equations (17), (18), (20)) as a function of the mass of 
the CP-odd Higgs boson, Ma- The other parameters used are those given in Eq. (40). 6 
As can be seen, the effect of applying propagator corrections at loop level as well as 
at leading order is significant. To obtain reliable results for the effective couplings the 
propagator corrections should therefore also be included at loop level. The curves for 
Options 1-3, where different methods are used to eliminate the IR divergences, are 
very similar to one another - i.e. there is little practical difference between the three 
approaches. 

In the following, we choose to use Option 3 - we make use of an altered coupling between 
the Higgs boson and a pair of charged Goldstone bosons and add the real radiation given 



5 Note that this only involves the CP-eigenstate Higgs boson H, as we are treating this particular 
diagram strictly at one-loop level, with no Z-factors applied. 

6 Plots varied over different parameters, and around a different point in the MSSM parameter space, 
show very similar features. 
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Figure 2: Numerical comparison of Options 1, 2 and 3 for dealing with IR divergent diagrams. 
Shown are the corrections to the form factor A for the WW final state as a function of Ma- 



by Eq. (10) (with the appropriate MSSM coupling constants) to render the process IR- 
finite. By using this option, we can calculate the ratio of NLO form factors 

1 . 2 Re [A* mp B (AA virt +AA 7 )] 
lAmp.NLol _ |Amp.B| + |A imp . B | 2 + ° real ^ 



l^ol 2 lA^eP , , 2Re[(AfM)*(AA^ t+ AAHM)] 

r T iSM 12 T" O; 



'real 



where AA 7 is the correction to the form factor from the virtual diagrams involving a 
photon (shown in Figures 1(a), 1(b)) and A/L v irt symbolizes all other virtual diagrams, 
so that 

AA, + AA virt = AAmp • (22) 

Using the relations between the SM and MSSM form factors, we can see that, when the 
effective HGG coupling is used: 

Amp.B = -^MSSM^tree 

A A 7 = FmssmAA^ 

2 Re jy!* mp B AAy] _ 2Re [(A™)' AA™] 

14- D |2 ~~ I4 SM |2 - W 

I^Mmp.BI I^Hreel 

Expanding, the ratio becomes 

|Amp.NLo| 2 |Amp.B| 2 / 1 + 5 7 + 5 v irt + <5rcal 



l^ol 2 IA S r ^l 2 \l + S-r + S^ + 5 Teal 



|2 



+ , (24) 

I ^tree I 
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where 

2Re[^ mp , B AA virt ] 2R^4^AA|a 

Ovirt — I 4 12 ' "virt— \ASM\2 " 

l^imp.BI l^treel 

Note that contributions from IR divergent diagrams cancel in the ratio of form factors. 
This cancellation happens only when we use Option 3, i.e. when we use an effective 
HGG coupling to eliminate two- loop IR divergent terms 7 . 

To include the contribution of the formfactor B we simply note that the ratio of polar- 
isation sums of the AB interference term and the A 2 term (cf. Eq. (6)) is 



/3 2 (l+lo 2 ) 
3/3 4 -2/3 2 + 3- 



r WW — oo2 i o m H l 2t) J 



with f3o from Eq.(ll). There is no interference between the formfactor C and the other 
two form factors. The ratio of partial widths may thus be written as 

r \A- B | 2 

P WW = FSM = I ASM 12 i 1 + ^ virt ~~ ^virt + r WW (Kirt ~ Kirt )] ( 27 ) 
1 I ^tree I 

with 

, 2Re[A* mp , B B] 2Re[(^ e )^ SM ] 

°virt i a i 2 ' °virt \ ASM \2 ' V > 

l^-imp.BI l^treel 

An analogous relation holds for the ratio of H — > ZZ partial widths. We have of course 
checked that the full ratio, including all contributions from photon diagrams, agrees 
well with the ratio as given in Eq. (24), and that all divergences (IR and UV) cancel in 
the full expression for the partial width. 

For off-shell decays the ratio pvff (V — W, Z) can be calculated with only a few 
modifications to the expressions above. First of all, the form factors A and B have 
to be calculated with one external gauge boson mass replaced by M/y/ (the invariant 
mass of the fermion pair). Furthermore, the ratio r V ff of polarisation sums has a more 
complicated form: 

= {Ml - {M v + M fr f){Ml - {My - M ffl f ){Ml — My — M] fl ) 
rvff 2(Mf I + M v + 10M v Mj f , + Mf fl -2Mf I (M v + M 2 fff )) ' [ ' 

With these modifications, the expression Eq. (27) can also be used for off-shell de- 
cays. In particular, the contributions from IR divergent diagrams still cancel in the 
ratio of (differential) partial widths, as long as Option 3 is used for the combination of 
propagator-type corrections and vertex corrections. 



7 Of course the mass of the SM Higgs boson My SM must be set equal to the heavy MSSM Higgs 
boson mass Mu- 
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3 The Parameter Scan 



As previously stated, important higher-order corrections to H — > WW, ZZ decays come 
from self-energy corrections of the initial-state Higgs boson (i.e. from the Z-factor con- 
tributions). These self-energy corrections also modify the Higgs boson mass and are 
dominated by loop diagrams involving the top Yukawa coupling, i.e. loops of top quarks 
and squarks. Beyond leading order, the H — > WW, ZZ decay rates therefore depend 
mainly on those MSSM parameters that enter the Higgs-stop-stop couplings and the 
stop mass matrix. These parameters are tan (3, the Higgsino mass parameter /i, A t (the 
stop trilinear coupling), mr (the right-handed stop mass) and m ? - i3 (the soft mass of 
the left-handed third-generation squarks). 

For small tan f3 the production of MSSM Higgs bosons at the LHC proceeds mainly 
through the loop-induced process gg — > H. In the SM this process is mediated by a 
quark loop. In the MSSM squarks can also appear in the loop. Again, the relevant 
contributions come from the diagrams involving the top Yukawa coupling. For large 
tan (3 the coupling of H to down-type quarks is enhanced by a factor 1/ cos (3, which 
can make the bb — >■ H production mode dominant [34-36]. In any case, the relevant 
MSSM parameters for the most important Higgs production processes are also tan/3, \i, 
A t , m~ tR and rriq L3 . 

These parameters are constrained by a number of experimental bounds. The strongest 
constraint comes from the observation of a Higgs-like resonance [1, 2], which we take 
to be the lightest MSSM Higgs boson. As the MSSM tree-level relations predict a light 
Higgs mass below the Z mass, the loop corrections have to push this mass up to 126 GeV. 
The dominant contributions to the Higgs mass still come from top and stop loops, but 
contributions from other sectors can also be relevant for satisfying the experimental 
bounds. Hence, the observation of a light Higgs boson at 126 GeV constrains many 
MSSM parameters simultaneously in a non-trivial way. Thus, to include the bounds 
correctly, we also have to consider MSSM parameters that have no significant impact 
on the heavy Higgs production or decay rates. These parameters include the gaugino 
masses and the soft masses of superpartners of the light fermions. Maximising the pp — > 
H — > WW, ZZ cross sections in the allowed part of this high-dimensional parameter 
space "by hand" would be both difficult and error-prone. We therefore rely on the 
numerical method of adaptive parameter scans as suggested in [37]. In the following 
paragraphs we describe our setup and the scanning method in detail. 

To determine the largest possible pp — >■ H — >■ WW, ZZ cross sections within the MSSM 
we scan over the following set of independent parameters: 

• tan (3, fi and Ma, 

• the gaugino masses M 1: M 2 and M 3 , 
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• the stop trilinear coupling A t 8 , 

• a universal soft mass mj for sleptons and sneutrinos, 

• a common soft mass rriq for all squarks except for the right-handed stop mass m^, 

• the mass of the light top squark, obtained by a judicious choice of the soft 
mass m ( - fl for right-handed stops after all other parameters are fixed. 

This scenario is a subclass of the so-called phenomenological MSSM [38]. It assumes all 
parameters of the soft SUSY Lagrangian to be flavour-diagonal, which is justified by the 
fact that new flavour structures in the MSSM Lagrangian are strongly constrained by 
flavour physics and thus have no significant effect in Higgs physics. Searches for SUSY 
particles at the LHC indicate that the superpartners of light quarks must be heavier 
than approximately 1 TeV. For simplicity, we use a common mass scale rriq for the 
corresponding soft masses. The LHC bounds on slepton and sneutrino masses are much 
weaker, so we use a different scale, mj, for the soft masses of sleptons and sneutrinos. In 
fact, if the anomalous magnetic moment of the muon — is used as an experimental 
bound, it is crucial to have a lower mass scale for sleptons and sneutrinos. Furthermore, 
the current LHC data cannot exclude a top squark that is lighter than the top quark. 
Thus, we keep the parameters A t and mi R independent, since they only enter the stop 
mass matrix. For convenience, we then trade the soft mass mi R of the right-handed stop 
for the physical mass of the lightest stop. 

Important constraints on this 10-dimensional parameter space come not only from direct 
Higgs searches at LEP, Tevatron and LHC, but also from the anomalous magnetic 
moment a M of the muon, the branching ratio BR(£> — > X s ^) and electroweak precision 
observables such as the p parameter, the effective leptonic mixing angle 9f s and the 
W mass. Thus, we discard portions of the parameter space according to the following 
criteria: 

• We discard any set of parameters that is excluded by direct Higgs searches at LEP, 
Tevatron and LHC at 95% CL. To do this, the masses and decay widths of the 
MSSM Higgs bosons and their effective couplings to SM fermions and gauge bosons 
are calculated with FeynHiggs 2.7.4 [12-18]. The effective couplings are (for the 
most part) implemented in the improved Born approximation, i.e. using higher- 
order corrections to the Higgs self-energies but no vertex corrections. 9 The masses, 
decay widths and effective couplings are then passed to HiggsBounds 3.8.0 [39, 
40], which confronts this information with 426 different search channels at LEP, 
Tevatron and LHC [41-95]. Internally, HiggsBounds uses a number of Standard 
Model results for the Higgs sector [96-128] to convert between experimental limits 
with different normalisations. Parameter sets that are excluded at 95% CL by any 
of these searches are discarded. 

8 The trilinear couplings At and A T are set to A t , and the trilinear couplings for the first and second 
fermion generations are set to zero 

9 Details about this approximation can be found in [15]. 
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To account for the discovery of a resonance at 126 GeV [1, 2] we discard any 
parameter sets for which the light MSSM Higgs boson mass lies outside the interval 
between 123 GeV and 129 GeV. 

As pointed out earlier, searches for supersymmetric particles at the LHC have 
already put strong constraints on the masses of squarks and gluinos. However, 
the interpretation of the individual searches in the context of a generic MSSM 
scenario is far from trivial and beyond the scope of this work. Thus, we simply 
require the gluinos and squarks of the first two generations to be heavier than 
1 TeV. Note that the squark mass limits do not apply to light stops due to their 
different production mechanism and decay pattern. For m x o > 80 GeV there 
is currently no lower limit on from LHC. For top squarks and all uncoloured 
supersymmetric particles we therefore use (largely model-independent) mass limits 
from LEP and Tevatron, as detailed in the SUSY review of [129]. Specifically, we 
require that 

- all slepton masses are larger than 100 GeV, 

— all chargino masses must be larger than 90 GeV, 

— the gluino and all squark masses except must be larger than 1 TeV, 

- m h > 100 GeV. 

The MSSM contributions Aa M to the anomalous magnetic moment a M = (g — 2)^/2 
of the muon are compared with the discrepancy between the SM prediction and 
the experimental value. We require [129] 

2 x 10~ 10 < Aa M < 36 x 10~ 10 . (30) 

The 2a range was extended by the uncertainty of the SM prediction for a^. 
The MSSM contributions Ap to the p parameter are restricted to [129] 

-0.0007 < Ap < 0.0033 . (31) 

After including MSSM corrections, the effective leptonic mixing angle 9f s is re- 
quired to satisfy [129] 

0.2280 < sin 2 6f < 0.2352 . (32) 

With MSSM corrections included, the W mass is required to satisfy [129] 

80.358 GeV < M w < 80.482 GeV . (33) 

For the branching ratio BK(B — > X s , -f) we impose very conservative limits. Our 
reasoning is that, unlike the other low-energy observables, BR(B — > X s ^) is very 
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sensitive to small violations of the assumption of minimal flavour violation. By- 
introducing new flavour structures in the MSSM we can therefore compensate 
deviations of BK(B — > X s 7) from its experimental value while leaving all other 
observables essentially unaltered. Thus, we conservatively require 

2.5 x KT 4 < BR(6 ->■ si) < 5.5 x 1(T 4 . (34) 

The correlation between the pp — > H — > VV cross sections and BR(6 — > 57) will 
be discussed in more detail later on. 

• In the vicinity of two-particle thresholds fixed-order calculations are numerically 
unstable and resummation techniques are required to obtain reliable results. To 
avoid this kind of numerical instability we discard all parameter sets where the 
mass of the Higgs boson H is within 2 GeV of the sum of masses of two particles 
that it couples to directly. This is not a physical constraint but merely a precaution 
to stop the scan algorithm from running into regions where our calculations are 
not reliable. 

MSSM corrections to the low-energy observables are calculated with FeynHiggs 2.7.4. 

As mentioned above, we use the numerical method described in [37] to systematically 
search for regions of the (experimentally allowed) parameter space where Ryy (as de- 
fined in Eq. (3)) is large. To this end, we must define an importance junction G(x), 
which is a real-valued function of the unknown model parameters: 

x = (tan f3,M A ,M 1 ,M 2 ,M 3 , A t , //,m,-, m^) . (35) 

We then use the VEGAS algorithm [130] to compute the integral of G over x. For the 
numerical integration we employed a modified version of the OmniComp-Dvegas package 
[131], which facilitates parallelised adaptive Monte Carlo integration and was developed 
in the context of [132, 133]. With each call to the integrand function, the parameters x, 
the value of Ryy and other relevant quantities are written to a file. This data may then 
be used to study the allowed range of Ryy and its correlation with other parameters and 
observables. If the importance function is chosen in a suitable way, the adaptive nature 
of the VEGAS algorithm guarantees that the "interesting" regions of the parameter 
space, i.e. those which exhibit relatively large values of Ryy, are sampled with a higher 
density. Note that neither the function G nor its integral have any physical meaning. 
The sole purpose of the importance function is to drive the adaptation of the VEGAS 
algorithm into those regions of parameter space we are interested in. An obvious choice 
for G(x) would therefore be 

{Ryy(x), x satisfies constraints 
0, otherwise 

where the 'constraints' are those discussed earlier in this section. The performance of 
the algorithm can be improved by using exponential dampening instead of "hard cuts" 
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m H [GeV] 

Figure 3: Scatterplot of Ftyy B (Eq. (3) evaluated in the improved Born approximation) 
against the heavy scalar Higgs mass Ma- 
in the importance function. Thus, we write G{x) as 

G{x) = RvvWllviiOiix)) , (37) 

i 

where the Oi denote all the observables used to constrain the parameter space. The 
functions yi are chosen to be equal to one for allowed values of the corresponding ob- 
servable Oi and to drop off exponentially outside the allowed range. In our scan, the 
constraints from the low-energy observables Aa M , Ap, sin 2 9f s , M\y and BR(6 — y 57) 
are implemented in this way. The constraints from Higgs searches are also treated in 
this manner. The related 'observable' is the ratio £95 of the signal cross section di- 
vided by the observed 95% CL limit for the most sensitive search channel, as provided 
by HiggsBounds. The mass bounds on SUSY particles, on the other hand, are imple- 
mented as hard cuts. Thus, our importance function may be non-zero for some points 
that do not satisfy the constraints discussed earlier. However, the scatter plots we show 
in this paper only contain points that satisfy the constraints. 

4 Results and discussion 

In the first part of this section we discuss scatterplots of Ryy B , as defined in Eq. (3) and 
evaluated in the improved Born approximation, against different input parameters and 
observables. The density of the points in these plots has no statistical interpretation. 
It is, however, safe to say that regions with a very low density of points can only be 
realised with rather finely tuned parameters. 

Figure 3 shows a scatterplot of Ry"^' 3 against the heavy scalar Higgs mass M#. We see 
that typical values of R l yy B do not exceed 0.15. On the other hand, the experimental 
results for SM H — y WW and H — y ZZ searches presented at the ICHEP2012 conference 
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Figure 4: Scatterplots of Ry"y' (Eq. (3) evaluated in the improved Born approximation) 
against the parameters A t , tan/3, mr and the quantity X t = A t — /xtan/3. 



probe, for certain Higgs masses, values of -R 1 ^' as low as 0.1. 10 In other words, the 
searches for heavy Higgs bosons already rule out certain regions of the MSSM parameter 
space. For M H < 160 GeV the values of Ryv drop rapidly below 0.05. The largest values 
of Ry^- B are reached for M H between 160 and 200 GeV. For M H > 240 GeV the typical 
size of R l yy B drops back below 0.1. 

In Figure 4 we show scatterplots of R l yy B against different input parameters. The 
quantity X t is defined as 

X t = A t - fxtan/3 , (38) 

so that m t X t is the off-diagonal element in the stop mass-matrix. We see that values 
of R l yy B above 0.1 are only possible for A t > 2TeV, tan/3 between approximately 
5 and 10 and X t > 2.5 TeV. Values of Ry^y' 3 larger than 0.15 additionally require 
m fl <150GeV. 

From Eq. (3) we know that Rvv is proportional to three ratios: the ratio cr^ SSM / erf^ of 
production cross sections, the ratio /r^ SSM of total decay widths (note the reversed 

10 Note that the constraints from direct Higgs searches implemented in our scan (via HiggsBounds) 
are only based on 2011 Higgs data. 
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order), and the ratio pyy of partial H — > VV decay widths. The last two quantities 
are universal in the sense that they appear in the signal strengths for all final states. 
Scatterplots of these ratios and the product of the two universal ratios against the value 
of R l yy B are shown in Fig. 5. We see that for values of Ry°y' above 0.1 the value of pyy 
usually does not exceed 0.05 and the remaining enhancement comes from the product 
of universal factors, which, for R l yy' B > 0.1, typically lies between 2 and 4. This is an 
interesting feature, since this enhancement factor would also appear in other decays like 
H — > tt or H — )■ bb. 

Fig. 6 shows the correlation of Ry^y' 3 with the two most constraining low-energy observ- 
ables, BR(6 — » sj) and Aa^. We see that large values of Ry"^' 3 typically coincide with 
large values of BR(6 — > sj) and small values of Aa M . The branching ratio BR(6 — > s'j) 
can be driven to smaller values by the chargino-stop loop contribution. This happens if 
the product pA t is small and negative. On the other hand, large values of R X yy B require 
A t to be large and positive and \p\ is bounded from below through the chargino mass 
limit. Thus small values of BR(6 — > sj) do not coincide with large values of R X yy B . 
A similar argument can be given for the correlation between -R 1 ^^ and Aa M . The 
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value of Aa^ is approximately proportional to tan/3/m?. Since large values of R^yy B 
require small values of tan and mj is bounded from below by the slepton mass limit 
we do not find any points where R^yy^ and Aa^ are large simultaneously. It should be 
stressed that, while these arguments make the features in Fig. 6 plausible, they do not 
suffice for a quantitative explanation. The correlations in Fig. 6 are really a result of 
the combination of different experimental bounds. 



To quantify the importance of the vertex corrections, we define the K factors Kyy as 

^Vr = ^riT • (39) 



imp.B 

Pvv 



where p l yy B denotes the ratio pyy from Eq. (3) in the improved Born approximation 
(i.e. with Z factors included but without any vertex corrections) and p v \ is the fully 
corrected value. Fig. 7 shows scatter plots of Kww and Kzz against the ratio Ry°y' B , 
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Figure 7: Scatterplots of the K factors Kyy (V = W, Z) due to genuine vertex corrections 
against the ratio i?™^' B . 
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m H [GeV] m„ [GeV] 

Figure 8: Scatterplots of the K factors Kyy (V = W, Z) due to genuine vertex corrections 
against the physical Higgs mass M#. Only points with Ry^y^ > 0.05 are included and points 
with iC/ > 0.1 are shown in black. 



i.e. the ratio Ryy calculated in the improved Born approximation. We see that, for 
Ryy ' B > 0.1 the K factors due to genuine vertex corrections typically lie between 0.8 and 
1.0. For R^ B > 0.05 the K factors can lie between 0.7 and 1.1. For R™* M < 0.02 the 
vertex corrections are numerically as important as the results obtained in the improved 
Born approximation. This means that, once the experimental sensitivity reaches two 
percent of the SM signal, the vertex corrections have to be included to derive limits on 
the MSSM parameters. 

Fig. 8 shows scatterplots of K ww and K zz against the physical Higgs mass M H . Only 
points with R vv > 0.05 are included and points with R l yy B > 0.1 are shown in black. 
We see that vertex corrections with a magnitude of more than 10% of the result in the 
improved Born approximation appear for Mh > 200 GeV and are always negative. The 
spike near Mh = 180 GeV originates from the ZZ threshold. 

The vertex diagrams contributing to H — > WW and H — > ZZ decays can be split into 
three separately IR and UV finite subsets: 

fermion/sfermion contributions from loops containing SM fermions or their super- 
partners, 

chargino/neutralino contributions from loops involving charginos or neutralinos and 

two-Higgs-doublet contributions from loops containing only gauge, Goldstone and 
Higgs bosons. 

To study the different contributions in more detail let us select a representative set of 
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parameters with a relatively large value of Ryv- 



m~ l = 300 GeV 
fi = 



M A = 260GeV , tan/3 = 7.5 , A t 
rrig = 2000 GeV , m i± 
2500 GeV , M 1 = -100 GeV 



150 GeV , m~ n 



1700 GeV , 

o = 1200 GeV 
M 2 = 200 GeV 
R zz = 0.0813 . 



=4> M^ = 246.7GeV , R ww = 0.0854 , = 0.0813 . (40) 

The vertex corrections to an d -Rzz are shown in Fig. 9 as functions of M A and 




200 
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M A [GeV; 



Figure 9: Genuine vertex corrections to Ryy (V = W,Z) as functions of M A (left column) 
and tan/3 (right column). All other parameters are fixed according to Eq. (40). The plots in 
the upper row show the ratios Rww ( re d lines) and Rzz (blue lines) with vertex corrections 
included. The dashed black line shows the same ratios in the improved Born approximation. 
{Rww an d Rzz are equal in this approximation.) The plots in the middle row show the K 
factors due to fermion/sfermion vertex diagrams (dashed lines), chargino/neutralino diagrams 
(dashdotted lines) and two-Higgs-doublet model diagrams (dotted lines). The solid lines com- 
bine all three types of corrections. Again, red and blue lines represent the WW and ZZ final 
states, respectively. The dependence of the mass Mh on Ma and tan (3 is shown in the plots 
at the bottom. The lines are interrupted if the corresponding parameter point violates the 
constraints discussed in Sec. 3. 



tan/3, with all other parameters fixed according to (40). We see that Rww an d Rzz 
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Figure 10: The off-shell K factor K z ^j as a function of the invariant mass M^j of the fermion 
anti- fermion pair // for the parameters of Eq. (42). The corresponding Higgs mass Mh is 
170.4 GeV, which is below the ZZ threshold. 



are different after the inclusion of vertex corrections, but the K factors only differ by 
a few percent. The dominant vertex contributions come from the fermion/sfermion 
diagrams, in particular from the diagrams involving top quarks and stops. For the 
chosen scenario these contributions to the K factors lie between —30% and —40%. The 
vertex corrections from the two-Higgs-doublet and chargino/neutralino sector typically 
only amount to a few percent. For Ma ps 260 GeV the H — )• hh threshold is crossed, 
which leads to the characteristic kinks in the graph. All three types of vertex corrections 
stem almost entirely from corrections to the form factor A (see Eq. (6)). Contributions 
from the loop- induced form factor B to Rvv are of the order of 10~ 4 and negligible for 
all practical purposes. 

Let us now examine the validity of the effective coupling approximation in the case of off- 
shell decays. As explained in Sec. 2.3, we can define the off-shell K factors K V ff(Mff/) 
by 

Pvff'(M f f>) 
P7Fr(M ff/ ) 

where V = W, Z, ff may be any SM fermion pair into which V can decay and Mfp is 
the invariant mass of the //' pair. If Ky/f is independent of Mffi (to a good approx- 
imation) the integrated MSSM cross sections for pp — >■ H — >■ V ff can be evaluated 
by scaling the corresponding SM cross section with Ryff (cf. Eq. 5) evaluated at some 
fixed value of Mff>. Such an approximation could be useful for Mh below 2Mz- To 
study the quality of the approximation we choose a scenario where M H is below the ZZ 
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Figure 11: Dependence of the ratio pw OI partial widths on the common complex phase 
4> of the parameters M\, M2, M3, p and At. Absolute values and relative signs are chosen 
in such a way that we re-obtain the parameters of Eq. (40) for = 0. The solid lines show 
the ratios pww (red) and pzz (blue) with vertex corrections included. The dashed black line 
shows the same ratios in the improved Born approximation. 

threshold: 

M A = 170GeV , tan/3 = 8.5 , A t = 2500 GeV , 
m t = 250 GeV , m~ q = 2200 GeV , m h = 250 GeV , m~ g = 1000 GeV , 
p = -900 GeV , Mi = -200 GeV , M 2 = 500 GeV 

M H = 172.2GeV , R ww = 0.102 , R zz = 0.100 . (42) 

Fig. 10 shows K Z fj as a function of Mjj. We see that K Z jj is constant within one per 
mille over a large range of Mffi. Scaling integrated SM cross sections for pp — > H — > 
Vff with Rvjf using the approximation of Eq. (5) therefore gives the correct MSSM 
cross sections with a relative accuracy of approximately one per mille. 

Finally, we would like to make a few remarks regarding the complex MSSM. As men- 
tioned earlier, introducing CP- violating phases in the MSSM Lagrangian leads to mixing 
between the CP-even and the CP-odd neutral Higgs bosons. We have seen that there 
is no interference between the form factor C and the (tree-level) form factor A. Thus, 
the largest values of pw are obtained if the decaying mass eigenstate has no CP-odd 
component, i.e. in the case of vanishing complex phases. This effect can be seen in 
Figure 11, where we take the parameter values of Eq. (40) and introduce a common 
complex phase </> for Mi, M2, M3, /1 and A t . Absolute values and relative signs are 
chosen in such a way that we re-obtain the parameters Eq. (40) for = 0. The plot 
shows the dependence of the ratio pw of partial widths (V = W, Z) as a function of (p. 
The dependence on can be much weaker for different choices of MSSM parameters, 
which only admit smaller mixing between the CP-even and CP-odd states. We have 
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run the full parameter scans with the complex phase as an additional variable, but 
the largest values of pvv were achieved for <p = 0. 

5 Conclusion 

We have calculated the partial decay widths for the heavy scalar MSSM Higgs boson H 
decaying into WW and ZZ final states at one-loop order and confirmed the available 
results of [19], which included the full one-loop corrections to the H — > ZZ decay 
width and the (s)fermion corrections to the H — > WW decay mode. We have extended 
the calculation to include the full 1-loop corrections in the H — > WW channel. To 
improve the precision of the one-loop result we proposed a method for combining Higgs 
propagator-type corrections (as defined in [15, 24] and calculated by FeynHiggs) with 
the genuine full one-loop vertex corrections for both H — > ZZ and H — > WW. We 
addressed the issue of infrared divergences appearing in the H — > WW process and 
ensured that our method leads to an IR finite result. In particular, no IR divergent 
diagrams need to be evaluated in the computation of the MSSM/SM ratio pvv (V = 
W, Z) of partial H — > VV decay widths if we use a modified coupling between the Higgs 
boson and a pair of charged Goldstone bosons. The same method allows us to calculate 
the MSSM/SM ratios pvff{Mffr) of differential partial widths for single off-shell decays 
H — > VV* — > Vff, where / and /' are two massless SM fermions and Mff/ is their 
invariant mass. We find that pvff is independent of M/y/ with a relative accuracy of 
approximately two per mille. Partial widths for single off-shell decay in the MSSM can 
therefore be safely estimated by scaling the corresponding (off-shell) SM partial widths. 

The possible size of the MSSM/SM ratios Ryv of Higgs production cross sections times 
branching ratios have been studied in an adaptive parameter scan. Experimental con- 
straints from several low-energy observables and direct Higgs searches at LEP, Tevatron 
and LHC were included in our scan (with the help of FeynHiggs and HiggsBounds). 
No assumptions about the SUSY breaking mechanism were made. We find that Rvv 
ratios of up to 0.2 can still be compatible with experimental constraints from direct 
SUSY searches and low-energy observables for M H > 160 GeV. These parameter space 
regions are currently or will soon be probed by the direct Higgs searches at the LHC. 
The one-loop vertex contributions to the decay processes typically lead to corrections 
between —30% and +10% for MSSM parameters where Rvv is larger than 0.05. For 
Rvv ^ 0.02 the vertex corrections can be numerically as important as the tree-level re- 
sults and Higgs self-energy corrections and therefore have to be considered when setting 
limits on the MSSM parameter space. 

The source code for our calculations is available on request from MW. 
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